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1 . Introduct ion . Definite evidence for particle acceleration in the 
solar wind came around a decade ago (Me Donald et al , 1976, Barnes et al, 
1976). Since then a considerable amount of data have been taken at dif- 
ferent distances from the sun confirming and extending the first 
findings and it is now widely accepted that particles are accelerated in 
the Solar Wind (Quenby, 1 983 ) . Two likely sources are known to exist: 
particles may be accelerated by the turbulence resulting from the super- 
position of AlfvSn and Magnetosonic waves (Statistical Acceleration) or 
they may be accelerated directly at shock fronts formed by the interac- 
tion of fast and slow solar wind (CIR's) or by traveling shocks due to 
sporadic coronal mass ejections. Naturally both mechanisms may be ope- 
rat ive. 


Shock acceleration has been widely investigated theoretically 
(Axford, 1981 and references therein) and there is substantial evidence 
of this as an operative mechanism in the Heliosphere (Scholer, 1984). 
However there are also experimental observations not obviously consistent 
with shock acceleration v.g.(Van Ness et al , 1984) and some ion enhan- 
cements only explained in terms of statistical acceleration (Richardson. 
1984). 


Previous treatments of statistical acceleration involved 
theoretical methods based on power spectral representations of the elec- 
tromagnetic field (see Quenby, 1983 and references therein). 

In this work the acceleration problem was tackled numerically 
using Helios 1 and 2 data to create a realistic representation of the 
Heliospheric plasma as will be described in the next section. Two 24 
hour samples were used: one where there are only wave like fluctuations 
of the field (Day 90 Helios 1) and another with a shock present in it 
(Day 92 of Helios 2) both in 1 976 during the STIP II interval. Transport 
coefficients in energy space have been calculated for particles injected 
in each sample and the effect of the shock studied in detail. 

2. Interplanetary Medium Model. The magnetic field used is defined from 
Bj three dimensional vectors where each one corresponds to the 8 sec. 
measurements of Helios 1 or 2. Every point in space is also furnished 
with a solar wind velocity vector \/j where these are the corresponding 

velocities in the solar wind frame. Consequently because all parts of 
the solar wind are moving there is an electric field Ej—Vi x Bj 

associated to every point in soace. For every sample we have divided 
the space in a series of layers where both the electric and magnetic 

fields are constant. Trajectories are integrated based on analytical 
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solutions of the equations of motion for each layer. For more details 
see Moussas et a 1 , ( 1 982 ) * 

Because we require layers of roughly 1/20 of a cyclotron radius 
only 100 MeV protons have been used in this work, lower energies 
require smaller IMF sampling time to keep a reasonable number of layers 
per gyroradius. 

3. Calculation of Transport Coefficients in Energy Space. First and 
second order coefficients (D^. and D^) were calculated in energy space 

by two different techniques. One is based upon the construction of a 
steady state distribution by injecting particles with a single energy 
To (100 MeV in this case) and removing them when they reach any of the 
two preset boundaries "above" and "below" the injection point (see 
Moussas et al , 1982 for details). 

Another method is based upon the time evolution of the energy 
distribution. D_ is calculated via a least square fit of succesive 
distributions first order momenta vs. time. If we take the diffusion 
equation in energy, integrate for injection at T=T 0 , assume a first 
order Taylor expansion in energy for D^, and make a determination of 

the spread <(T t -<T>) 2 > we arrive at 

<(T-<T>) 2 > = 2D TT (T 0 )t+2(D T ) 2 


from where can be calculated. 

Both methods are generalizations of those initially developed 
by Jones et al (1978) to study pitch angle scattering on a randomly 
generated IMF. 


The effect of the shock is traced by counting every particle 
encounter with it and recording the resulting energy change. The 
average energy change is given by the ratio of the total gain and the 
number of shock encounters. Transport coefficients can be calculated 
using the time an average particle takes to get back to the shock 
t„ = 2A/v where A = mean free path = 0.03 AU (see Val d6s-Gal icia et al , 
1984) and v = particle velocity. 


4. Shock Characteristics. The shock used in this study passed 
through Helios 2 on day 92 of 1976 when the spacecraft was at 0.45 AU 


from the Sun. It is a perpendicular shock 
1971) with a high Alfv6n Mach number M^ = 7.5. 


0 Bft = ^9° ( Le PP' n 9 et a1 


The magnetic field 


overshoot is some 24 Gammas which is around 50 % of the downstream 
field. 


It was assumed to be a plane shock and data were transformed 
to a frame where the Y"Z plane coincides with that of the shock so 
that it is parallel to our layer planes. 
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5. Results and Discussion. In figure 1 we present a result of a time 
stationary distribution of particles for an experiment carried out 
using data from day 90/Helios 1. It can be easily appreciated the 
asymetry of the distribution towards the "right" of the injection 
energy T 0 showing strong acceleration by waves. Boundaries were put 
at + 2 MeV from T 0 - 100 MeV. It should be noted that the mean free 
path these particles is very small (A = 0.006 AU, Va1d6s-Gal icia et al , 
198A). 


Figure 2 shows an evolution diagram of the distribution 
function of particles in energy and time. Contours are drawn for 
different density levels. Areas with zero density within the distri- 
bution are shown black. Injection energy is at the middle of the 
figure and energy increases downwards. Full extension from top to 
bottom is 6 MeV. The horizontal scale covers 2000 sec. Balck diamonds 
represent first order momenta of particle distributions every 100 secs 

In table I we show the main results of this work. TransDort 
coefficients are average values of the calculations done by the two 
different methods. Although not shown individually both values agree 
within 15%. The second row of table I shows results for all particles 
( 85 ) used with Helios 2 day 92 sample including 79 shock crossings so 
they are representative of the two proccesses involved (shock and 
statistical acceleration). 

We can aprecciate that even though average energy changes 
corresponding to statistical and statistical plus shock acceleration 
are greater than the shock produced average, transport coefficients 
for this process are greater by an order of magnitude. In the last 
column of the table we have calculated the the e-folding times (t) 
corresponding to every experiment (Wibberenz et al, 1972). If we 
calculate the time for adiabatic deceleration at 0.^5 asuming 
V sw = 2,00 and radial expansion we get T gd = 35,3 hours. Thus 

al tough statistical acceleration can have an effect in reducing the 
adiabatic cool ing only the shock accelerated particles are able to 
overcome it at these energies. 

Unfortunately we were limited by the time resolution of the 
data and could not explore lower energies in a similar manner, 

6. Conclusions, 


a) The time evolution and time stationary methods to calculate 
transoort coefficients agree quite well. 

b) Statistical acceleration at particle energies of 100 MeV in the 
inner Heliosphere is not able to overcome the effects of 
adiabat ic cool ing. 

c) Shock acceleration as opposed to statistical acceleration (nay be 
a source of particles at these places an energies. 
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d) Statistical acceleration may be a more efficient process at 
lower energies. 




FIG I FIG 2 


TABLE I 


LOV/SC 

POSITION 

CAU) 

ACCELERATION 

PROCESS 

^aT> 

(r'EV) 

°T 

(MeVs" 1 ) 

d tt 

(He vs' 2 ) 

T 

(HR3) 

9C/HEl 1 

0.31 

etati:tical 

0.76 

(3.6+2.UX10' 4 

(4::3>xl0' 4 

76.5 

92/HEL 2 

0.45 

SH0CK+S7ATISTICAL 

0.83 

(5.2+2)xl0' 4 

C5+3,5)XlO' 4 

53. 

92/HEL 2 

0.45 

SHOCK 

0.26 

(3.4+1.5)xl0" 3 

(l+.9)xl0' 4 

10.9 
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